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H II regions are volumes of gas surrounding high-mass stars that are ionized 
and heated by the stellar ultraviolet radiation. This ionization and heating 
significantly raises the thermal pressure of the gas, which is the principal 
driver of the dynamics of all but the very smallest and very largest H II regions. 
The simplest model of H II region evolution predicts a slow and decelerating 
expansion of the ionized gas with very little internal motion, but this is rarely 
observed. Instead, when examined in detail, such regions are found to be 
highly structured with complex internal motions. 

For reasons of space, this review is rather narrowly focused on the tradi- 
tional H II regions that are found in our galaxy, ionized by one or a handful 
of O stars. The dynamics of photoionized gas has a much broader domain 
of application than this, covering such objects as planetary nebulae, novae, 
Wolf Rayet nebulae, broad and narrow line regions of active galaxies, and the 
reionization of the universe by the first generation of stars. Many of these are 
covered in other chapters in this volume. The physical priniciples presented 
here will still apply in these broader contexts, but care must be taken, as some 
of the shortcuts and approximations commonly used for H II regions may no 
longer be valid. 

The chapter is divided into three sections. The first section introduces 
the equations governing the dynamics and physical structure of H II regions 
and discusses the approximations that are commonly employed. The second 
section presents the broad physical concepts that provide the building blocks 
for constructing global models of Hll regions. The third section provides an 
overview of such models, as applied to the different evolutionary stages of 
Hll regions, with particular emphasis on recent models of the evolution in a 
clumpy, turbulent medium. 
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1 The equations 

The principal equations necessary for calculating H II region dynamics under 
the standard assumptions used in the field can be divided into three groups. 
First, the Euler equations, which describe the motion of the gas. Second, the 
radiative transfer equation, which describes the emission and absorption of 
photons. Third, the rate equations, which describe the transitions between 
different atomic/ionic/molecular states. For each group, I first present the 
general form of the equations, in which many difficult bits of physics are hidden 
away in innocuous looking terms, before discussing the sorts of approximations 
that might be made in different types of models (all attempts to theoretically 
model H II regions involve at least some level of approximation) . 



1.1 Euler equations 

These describe the motion of a non-viscous, non-relativistic gas of density p, 
pressure P, and vector velocit y u, a ll of w hich are functions of position, r, 
and time, t. See, for example, Shu (1992) for full derivations. The ratio of 
specific heats is assumed to be 7 = 5/3, as is appropriate for a monatomic 
gas. The equations are given here in there conservation form. Conservation of 
mass gives 

^ + V.(/m) = 0. (f) 
Conservation of momentum gives 

j t (pu) + V(P + pu 2 )=gp, (2) 

where g is the acceleration of the gas due to "body forces" caused by gravity, 
magnetic fields, or radiation pressure. Conservation of energy gives 
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u-gp + H-C, (3) 



where H and C are the volumetric rates of gas heating and cooling, respec- 
tively, due to atomic processes. These equations are supplemented by the 
ideal-gas equation of state, which defines the gas temperature, T: 

P=*°£, (4) 
pmu 

where fee is Boltzmann's constant, ijih is the mass of a hydrogen nucleus, and 
p is the dimensionless mean-mass-per-particle (for solar abundances, p ~ I .3 
for neutral gas and p ~ 0.6-0.7 for ionized gas, depending on the ionization 
state of helium) . 

Although the above single-fluid treatment is generally adequate for the 
gas, it is not always such a good approximation when dust grains are con- 
sidered in detail. In such a case, one can employ a multifluid treatment 
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Gail fc Sedlmavrlll979tM . in which one uses a separate momentum equa- 
tion for each fluid, including explicit interaction terms for collisions between 
particles of the different fluids. 



1.2 Radiative transfer 



The fundamental equation of radiative t ransfer describes the behavior of the 
specific intensity, I v {n,r) l)Mihalaslll978|) . which is a function of frequency, v, 
position, r, and direction, n: 

— rr + n ■ V I u = r\ v - Xulu, (5) 
c at 

where c is the speed of light, and n v , Xv are respectively the emissivity and 
absorption coefficient, which may contain an arbitrary amount of physics. For 
many purposes, it is sufficient to work with angle-averaged moments of the 
specific intensity, such as the mean intensity, 

47rJ„= f I v dQ, (6) 

J47T 

and the radiative flux, 

F v = [ fil v dQ. (7) 



1.3 Rate equations 

The general equation for the evolution of the partial density of a particular 
state, i, can be written as 

^ + V-imu) = Gi + J2 Rj^inj - nt I Si + ^ R^j , (8) 

where iZ,-_ >j, Rj~>i represent the rates of transitions between state i and other 
states j, whereas Gi, Si represent respectively sources and sinks of state i 
due to other processes. This form of the rate equation can apply to internal 
states of an atom, an ion, a molecule, or a dust grain. Equally, it may apply 
to the the total abundance of a given ion stage. In general, the rates may be 
functions of the local radiation field, electron density, temperature, or indeed 
any other physical variable. 



1.4 How to avoid the dynamics 

The most drastic simplification from the point of view of the dynamics is 
to assume that there is no dynamics at all. This is the static approxima- 
tion, which consists in setting u = and dX/dt = for all quantities X. 
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In this case, equation Q is trivially satisfied, while equation @ reduces to 
hydrostatic equilibrium, or P = constant if there are no external forces. The 
energy equation © reduces to H = C and the rate equations (JSJ) reduce 
to sources = sinks. This approximation is commonly employed in photoion- 
ization codes, which treat the microphysics of the energy and rate equations 
in great detail, and usually the radiative transfer also. As is shown below, 
this approximation can be an acceptable one for the ionization and thermal 
balance of the interior of an H II region in the many instances where the ion- 
ization/recombination and heating/cooling timescales are much shorter than 
the dynamic timescale. 

Slightly more realistic is the steady-state or stationary approximation, in 
which one allows for non-zero gas velocities, but still ignores all derivatives 
with respect to time. This adds the advective terms, of the form V-(tiX), to 
the static form of the equations. These terms represent the material flow of 
the quantity X from one point to another and are most important in regions 
of strong gradients in u or X. In many instances, the flow timescale through 
the part of the Hll region of interest is short compared with the timescale 
for changes in the parameters of the flow, due, for instance, to changes in the 
incident ionizing flux or the neutral density outside the ionization front. In 
such cases, the stationary approximation is often a very good one, provided 
that a suitable reference frame is chosen. In other cases, especially when the 
long-term, global evolution of the Hll region is of interest, or if dynamical in- 
stabilities are expected, then it is necessary to include the non-steady dX/dt 
terms and solve the fully time- dependent equations. Even in this case, it is usu- 
ally sufficient to treat the radiative transfer equation in the stationary limit, 
except whe n light-travel times are significant compared with other timescales 
of interest ( Shapiro et al.ll2005^ . 



1.5 How to avoid the atomic physics 

In other contexts, one may wish to study the dynamics in detail without wor- 
rying unduly about the finer points of the microphysics or radiative transfer. 
To that end, various rules of thumb have evolved, which give satisfactory 
approximations in many common situations. However, the use of these ap- 
proximations should be justified on a case-by-case basis. The simplifications 
presented here are at the extreme end of what one can get away with in a 
toy model. For more serious applications, one should consult a text such as 
lOsterbrock k, Ferlandl ( 2006^ for extra ingredients to add. 



Thermal balance 

The temperature in galactic H II regions is determined principally by the bal- 
ance between photoelectric heating (ejection of energetic electrons by pho- 
toionization) and cooling due to forbidden lines of the ions of heavy elements, 
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which are excited by electron collisions. The resultant equilibrium tempera- 
ture is ~ 9000 K and rarely varies by more than 50% throughout the region. 
It is therefore natural to use an isothermal approximation for the ionized gas. 

Strong dynamic effects can cause this approximation to break down, but 
this is rare in Hll regions. The dynamical terms in equation J2J will not 
be greatly important unless the dynamic timescale is shorter than the heat- 
ing/cooling timescale, which is typically 3-10 times shorter than the recombi- 
nation timescale. This requires very high gas velocities of > 500 km s^ 1 unless 
the ionization parameter (see below) is much lower than is typical found. 

Transfer of ionizing photons 

Since hydrogen is (usually) the most abundant element, it is natural to con- 
sider a hydrogen-only approximation, in which the only photons of interest 
are those capable of ionizing hydrogen from its ground state and the only 
opacity source considered for those photons is the photoionization process 
itself. If the diffuse field is treated in the on-the-spot approximation (see 
lOsterbrock k, Ferlandl I200II sec. 2.3), then only the direct stellar radiation 
need be explicitly considered. Assuming the presence of only one ionizing star, 
then the radiation field is monodirectional, so that all angular moments of the 
radiation field are equal in magnitude, and in particular 47rJ„ = F v = \F U \. 
In this case, one can get away with considering only the frequency-integrated 
flux of ionizing photons: / = J {F y jhv) dis, where vq is the frequency cor- 
responding to the Lyman limit. The radiative transfer equation then reduces 
to V-/ = — o~n n f, where r is the radial direction from the star, n n is the 
number density of neutral hydrogen atoms and a is the frequency-averaged 
photoionization cross-section, weighted by the local ionizing spectrum. In this 
approximation, a is a function only of r = J o~n n dr, and can be precomputed 
for a given stellar spect rum. In reality, dus t opacity will often m ake a signif- 
icant contribution (e.g.. lArthur et al.ll2004i lErcolano et al.ll2005h and should 
be included in any realistic treatment. 

Ionization balance 

Collisional ionization of hydrogen is unimportant for T < 20, 000 K and three- 
body recombinations are negligible at typical Hll region densities. The hydro- 
gen ionization balance is then simply between photoionization and radiative 
recombination, which in the approximation discussed in the previous section 
is 

— £ + V-(npu) = n n af - an p n c , (9) 

where n p , n c (~ n p ), and n n are the number densities of ionized hydrogen, 
electrons, and neutral hydrogen, respectively, and a is the appropriate recom- 
bination coefficient (Case B if the on-the-spot approximation is used). If the 
radiative transfer is treated in more detail, then af should be replaced by 
ffiu{^KJ v /hv) dv. 
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2 Physical concepts 

Before considering models for Hll regions as a whole, it is useful to break 
down the problem into distinct physical ingredients, which can be studied 
separately. 

2.1 Static photoionization equilibrium 

In order to compare and contrast with the results of dynamical models, 
it is instructive to first consider the static case. The algebra is simplest if 
one considers a uniform, plane-parallel slab, of density n, one face of which 
illuminated by a perpendicular ionizing flux, /o, and in which all the ex- 
treme simplifying assumptions of the previous sections are assumed to hold. 
The degree of ionization is defined as x — n p /n and, for reasonable val- 
ues of /o and n, this fraction is very close to unity at the illuminated face: 
1 — xq ~ an/(afo) ~ 3x 10~ 6 Y , where Y = /o/(cn) is the ionization param- 
eter, or ratio of photon density to particle density, which typically takes values 
from 10 -4 to 10 -2 . As one moves to greater depths, z, in the slab, then integra- 
tion of the ionization balance and radiative transfer equations shows that the 
ionizing flux decreases approximately linearly with depth, f(z) ~ fo(l — z/zo), 
until one reaches a depth Zq ~ fo/(an 2 ), where the ionization fraction swiftly 
falls from x ~ 1 to x ~ over a short distance, Sz ~ l/(n<r). This depth 
to the ionization front, Zo, is referred to as the Stromgren distance and it is 
apparent that Sz/z n ~ (1 — x ), confirming that the front is indeed thin. The 
ionization parameter, T, is the single most important parameter describing 
the H II region; as well as determining the degree of ionization and the relative 
thickness of the ionization front, it is also proportional to the column density 
through the ionized gas: nz = /o/(cm) ~ 10 23 Y cm~ 2 . 

This last result allows one to see how the ionized region in approximate 
photoionization equilibrium must respond to a change in density. If the density 
in the slab increases, perhaps as the result of compression by a shock, then 
the ionization parameter decreases and so must the column of ionized gas, 
meaning that some of the gas must recombine. Contrariwise, if the density 
decreases, perhaps due to expansion, 1 then the ionization parameter and the 
ionized column will increase, leading to the advance of the ionization front 
into previously neutral gas. 

2.2 Ionization front propagation 

The speed of propagation, U, of a plane, steady ionization front can be de- 
termined by considering the jump conditions in the physical variables across 



1 The same holds for expansion in a spherical geometry since n (~ r 3 ) falls more 
rapidly than / (~ r~ 2 ). 
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it teahnlll954l) . Variables on the neutral side of the front are given a sub- 
script '1', and those on the ionized side a subscript '2'. The velocity in the 
frame of reference of the ionization front of the neutral gas, V\, is equal in 
magnitude to U, but opposite in direction. It is apparent that one has /i ~ 
and (1 — 2:2) — 0, so that by integrating equation JSJ) across the front 2 one 
has vi = jiln\. When one then considers mass and momentum conservation 
across the front, assuming isothermal sound speeds of ai and a 2 on the two 
sides, one finds two classes of solution (for example. IShulll992T) . Fronts with 
v > v-r ~ 2d2 — 20 km s _1 move supersonically into the neutral gas and are 
called R-type (for "rare"). Fronts with v < vn ~ a 2 /(2a 2 ) ~ 0.05 km s _1 
move subsonically into the neutral gas and are called D-type (for "dense"). 
If conditions are such that f 2 /n\ falls between vr> and then the ioniza- 
tion front will be preceded by a shock, which compresses the neutral gas to a 
higher density, n[, so that / 2 /^i < «d and the ionization front becomes D- 
type. Fronts can be further divided into those that contain an internal sonic 
transition, which are termed strong, and those that do not, which are termed 
weak. The only stable R-type fronts are weak, in which the front moves su- 
personically with respect to both the neutral and the ionized gas and such 
fronts show a low density contrast between the two sides: 1 < ri\jn 2 < 2. 
The limit of extreme weak D-type fronts as v\ — > corresponds to the static 
case considered above, in which the gas pressure is constant across the front so 
that 711/712 = a\/a\ ~ 100. When vi = vu, the front is termed D-critical, with 
an ionized velocity that is exactly sonic with respect to the front, v 2 — a 2 , 
and an even higher density contrast, n\jn 2 — 2a\/a\ ~ 200. Strong D-type 
fronts can only occur if the sound speed has a maximum a m > a 2 inside the 
front and give a m < v 2 < 2a m . They are probably most relevant at very high 
metallicities, where the equilibrium temperature of fully ionized gas can be 
much lower than that of partially ionized gas. 

When a volume of gas is initially exposed to ionizing radiation, the flux 
is usually high enough that the ionization-front is R-type and propagates 
supersonically through the gas. However, as the front progresses to greater 
distances, an increasing proportion of the flux is used up in balancing recom- 
binations in the ionized gas. After a time of order l/(an) and a distance of 
roughly one Stromgren distance, f 2 has become low enough that U < 2a 2 and 
a preceding shock detaches from the ionization front, which becomes D-type. 
Soon after this point, the speed of the ionization/shock front falls below the 
sound speed in the ionized gas, so that the subsequent evolution of the front 
becomes sensitive to the internal conditions in the H II region, which deter- 
mine the boundary conditions on the ionized side of the front. In particular, 
if the ionized gas is free to flow away from the front, then the front is liable to 
remain approximately D-critical. On the other hand, if the ionized gas is con- 
fined, perhaps by a closed geometry (nowhere for the gas to go to) or a high 
imposed pressure (e.g., from a stellar wind), then the front will be weak D. 



2 In this approximation, recombinations in the front are ignored. 
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The latter case corresponds to the classical evolution of a Stromgren sphere, 
which is described in many textbooks fe.g.. lDvson fc Williamdll997j) . 

The relative importance within the front itself of dynamic terms in the 
ionization balance is large, of order u^afc*. ~ 10M, where M — ^2/^2 is the 
Mach number of the ionized gas (M = 1 for a D-critical front). However, given 
that the front is generally thin compared with the H II region as a whole (sec- 
tion I2.l)l . the relative importance, A a d, of dynamics in the global ionization 
balance is much smaller, A a d = uo,/(Te) ~ 0.003-0.3. This is why photoion- 
ization models that ignore dynamics are often a reasonable approximation for 
the interior of an Hll region. 

2.3 Structure of a D-type front 

When examined in detail, an ionization front is found to have considerable 
internal structure. An example of an approximately steady D-critical front 
is shown in Fig. ^ which results from a radiation-hydrodynamical simulation 
that uses rather simplified atomic physics and an artificially lowered ionization 
parameter in order to resolve the ionization front on a two-dimensional grid 
ijHennev et al]l2005a|) . However, the principal results compare well with much 
more detailed one-dimension al calculations of wea k-D fronts using a state-of- 
the-art plasma physics code ijHennev et al ] l2005bl) . In order to emphasize the 
regions of the front where interesting changes occur, the physical variables are 
plotted as a function of the mean optical depth for ionizing radiation from 
the star. The velocity is in the direction of lower optical depths and is given 
relative to the ionization front, which is moving at 0.95 km s _1 away from the 
star. 



Net heating rate 
■■■■ Density, n/10 s cm" 3 

Pressure, P/10~ 7 dyne turr 2 

Temperature, T/10 4 K 

.... Velocity, vj km s -1 

Ionization, — x) 

m*m. Sound speed, cj km s -1 



0.1 0.2 0.5 1 2 5 10 20 50 100 

Optical depth to ionizing radiation, r 

Fig. 1. Physical variables in a typical D-critical ionization front taken from a cut 
along the axis of a two-dimensional hydrodynamical simulation of a champagne flow. 
Stellar photons are movi ng from left to right, but the gas is moving from right to 
left. Figure adapted from lHennev et, all l)2005al) . 
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The net rate at which the gas gains energy from atomic processes (pho- 
toelectric heating minus radiative cooling; solid gray line) has two separate 
peaks in the front. The deeper peak occurs at r ~ 8 and results in the heat- 
ing of the gas from 100 K to 9000 K and its acceleration from 0.1 km s _1 to 
3 km s _1 , but leaves the gas still largely neutral (x ~ 0.01). The shallower, 
broader peak occurs at r ~ 3 and results in the complete ionization of the gas 
and its acceleration to ~ 10 km s , with very little change in temperature. 3 
The gas then passes through a sonic point at r ~ 1, at which point x ~ 0.95, 
and continues accelerating gradually up to ~ 30 km s" 1 as it flows back past 
the star (not shown). Although the far downstream velocity of the ionized 
gas is supersonic, the front is still D-critical since there is no internal maxi- 
mum in the sound speed. Instead, it is the divergence of the gas streamlines, 
as the geometry deviates from plane-parallel, that allow the gas to continue 
accelerating once it is fully ionized and isothermal. 

The shock that precedes the ionization front is also visible, at depth of 
t ~ 50. The shock is propagating into the neutral gas slightly faster than 
the ionization front velocity, so this part of the model is not exactly time- 
steady, but rather the column of shocked neutral material grows with time. 
This part of the structure is the least well-modeled in the simulation, since 
the physics of the photo d issociatio n region is not treated adequately, but 
sec H osokawa fc Inutsukal i|2005bla|) for dynamical models that do treat this 
properly. 

2.4 Photoablation flows 

If the neutral density is similar in all directions from the ionizing star, then 
one would expect the radius to the ionization front to be also similar, which 
would lead to a global geometry for the front that is closed, or concave. Despite 
this, the opposite case is of great relevance to Hll region dynamics, in which 
the ionization front is convex; that is, it curves away from the ionizing star. 
This can arise, for example, if dense condensations exist in the neutral gas, 
which slow down the propagation of the ionization front, causing it to "wrap 
around" the obstacle. If this situation persists, then a steady flow of ionized gas 
will establish itself that accelerates away from the convex front. The density 
in this ionized photoablation flow increases sharply towards its base, where 
it may greatly exceed the mean density in the Hll region. Models for such 
flows were originally developed in the context of the p hotoevaporation of 
small, self-gravitating globules llDvsonlll968l lKahnl ll969) but have a much 
wider range of application. By adopting a cylindrical, rather than spherical 
geometry, similar models can be applied to cases where the ionization front 
wraps around a dense filament, as in the Bright Bar in the Orion Nebula, 



3 This can be compared with a simple analytic estimate for the optical depth to the 
ionization front (see Sec. 12.21 : ry = — lnA a d — 10.3 + 2.3 log 10 T, giving th ~ 1 
to 6 for ionization parameters in the range T = 10 -4 -10 -2 . 
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and they can even be extended to cases where the ionization front is flat, or 
slightly concave, as in the case of champagne flows (see section EOl below - ) . 
In all cases, the structure of the ionized flow is similar, provided only that 
the gas streamlines diverge and that any pressure imposed at the downstream 
boundary is low compared with the thermal pressure of the ionized gas at the 
ionization front. In order for a steady-state approach to be valid, it is also 
necessary that the timescale for significant changes in the configuration of the 
neutral gas be long compared with the dynamic timescale of the ionized flow. 
This condition is easily satisfied in some cases, such as the Orion proplyds, but 
only approximately satisfied in othe rs, such as the radiation-driven implosion 
of a neutral globule llBertoldilll98fll) . 

The simplest models proceed by patching together the solution for a plane- 
parallel D-critical ionization front (see previous sections) with an isothermal 
wind solution for the fully ionized flow. This is valid so long as the width of 
the ionization front is very small compared with its radius of curvature, r c , 
which is almost invariably the case. The gas leave s the ionization front at the 
sound speed but accelerates rapidly (see Fig. 6 of lHennevll200.'^ . so that the 
ionized density falls off more quickly than it would in a constant velocity wind, 
approximately mimicking an exponential profile. As the Mach number of the 
flow increases, the acceleration gradually lessens and a terminal velocity is 
reached when the isothermal approximation breaks down, which occurs when 
the expansion and heating timescales becoming comparable. For a spherical 
flow with A a d = 0.01, this occurs at velocity of ~ 40 km s _1 and a radius of 
~ 10r c . 

Unless the ionization parameter is rather low, A a( j will be small and the 
flows are recombination dominated, which means that the ionized flow has a 
high optical depth to ionizing radiation and is in approximate static ionization 
balance. 4 The effective recombination thickness of the flow 5 is h e g = 0.12r c 
for a spherical geometry or h e R = 0.3r c for a cylindrical geometry, so long 
as the curvature is positive (convex) and r c <C r+, where is the distance 
of the front from the ionizing star. For the case where r c > r*, one instead 
finds h c ff ~ 0.2-0.7r\t (see Fig. 7 of lHennev et al.ll2005a|h In either case, the 
entire flow structure is uniquely determined by the incident ionizing flux and 
the curvature of the front. If r c < 0.2r*, then the photoevaporation flow is 
not capable of evacuating the region all the way back to the ionizing star, so 
the incident flux may be reduced by recombinations in the interior of the H II 
region. 



4 The opposite adve ction dominat ed case is generally only seen in flows with low 
surface brightne ss jHejrngjtoO-l j) I , most notably the knots in the Helix planetary 
nebula llLopez- Martin et alJl20Q]llO'DeU e t alJ l200dh in which a high fraction of 
incident photons reach the ionization front to ionize new gas. 

5 Defined as n^h^s = J n 2 dr, where n is the ionized density with value no at the 
ionization front. 
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2.5 Other ingredients 
Stellar Winds 



High mass stars invariably drive fast (~ 1000 km s 1 ) stellar winds ( Lamers fc Cassinellil 



1999), Wind mass loss rates are still very uncertain du e to the poorly under - 
stood effects of clumping on observational diagnostics l|Fullerton et alJl2006|) . 
but are probably in the range 1O _9 -1O _7 M0 yr _1 for main-sequence O stars. 
The effect of a stellar wind on the dynamics of an H II region depends cru- 
cially on whether or not a hot bubble of shocked stellar wind can persist, 
which determines whether the shell of shocked Hll region gas is energy-driven 
or momentum-driven (see the chapter by Arthur chap ter in this volume). 
In the energy-driven case, ICapriotti fc Kozminskil l)200l|) show that the shell 
swept up by the stellar wind can dominate the expansion energy of the region 
at e arly evolutionary times and when the ambie nt density is > 10 4 cm -3 (see 
also lDvsor]ll977HGarcia-Segura fc Francolll996|) . The wind-blown shell is also 
capable of trapping the ionization front in some instances. In the momentum- 
driven case, the wind never dominates the global energetics of the region but 
can nevertheless have an important local influence on the internal dynamics. 
One example is the bowshocks foun d around proplyds in the inner Orion Neb- 
ula l|Oarcia-Arredondo et alJboOlh . Various mechanisms may prevent the for- 
mation of the hot shocked bubble t hat is required for the energy- driven case, 
including therma l conduction (e.g.. IPorland fc Montmerlelll987j) or dynam- 
ical instabilities ijBreitschwerdt fc KahrJll988|) at the contact discontinuity, 
or ma ss-loading of the wind due to embedded photoa blating clumps or pro- 



plyds IjDvson et al.lll995t iGarcia-Arredondo et alJl2002|) . All these can cause 
enhanced cooling in the shocked wind, leading to the loss of its thermal pres- 
sure. Specific instances of the importance of winds in the evolution of H II 
regions are described in section [3] below. 



Radiation pressure 

The transfer of momentum between the radiation field and the gas can some- 
times affect the dynamics of Hll regions. This arises in three different ways: 
first, through trapped resonance line radiation (mainly Lyman a); second, 
through the momentum of stellar photons transferred during the photoion- 
ization process; and, third, indirectly through collisional coupling with dust 
grains that are accelerated by the absorption of stellar radiation. 

Although the diffusion of resonance line photons can in principle intro- 
duce non-local couplings between disparate parts of an Hll region, this is 
sever ely hampered by the presence of dust absorption ((Hummer fc Kunaszl 
1980), which allows a purely local approach to be used if the medium is ho- 
mogeneous on scales < 0.1/^, where k <\ is the dust absorption coefficient. 
By this means, iHennev fc Arthur! l)l998|) showed that the resonance line ra- 
diation pressure is proportional to the gas pressure, making a contribution of 
only ~ 5% for a standard dust-gas-ratio. 
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The radiative acceleration due to absorption of ionizing photons of mean 
energy (hv) is g ra( \ ~ ((hv) /c)(f / ph^ s ), where ft. a b s is the thickness over 
which the photons are absorbed. The acceleration due to pressure gradients 
in a photoablation flow fsection l2.4|l is <?fl w — a? /h e g. Comparing the two, and 
putting /i a bs = h e s, one finds 3rad/ffflow — T(hv)/{2kT) ~ 10T, so that the 
pressure of the ionizing radiation is of only secondary importance unless the 
ionization parameter is higher than is typically found in H II regions. Assuming 
that the gas and grains are effectively coupled, one finds a similar result for 
the absorption of radiation by dust, although this can become more important 
for cooler stars, where the relative luminosity at non-ionizing wavelengths is 
higher. Radiation pressure on dust can also be important in the inner parts 
of H II regions due to the increased flux close to the io nizing source, 6 where 
it ma y be responsible for central cavities in some cases l)Mathewslll967t Ilnouel 
2002). 



Magnetic fields 

The magnetic field makes a significant, perhaps dominant, contribution to 
the pressure both in molecular clouds and in the diffuse ISM. Hll regions are 
generally over-pressured with respect to their undisturbed surroundings, but 
are within a factor of two of pressure equilibrium with the shocked neutral 
gas that surrounds them. The dynamic importance of magnetic fields inside 
the H II region depends sensitively on how the field strength, B, responds to 
compression in the shock and rarefaction at the ionization front. This can 
be approximately characterized by an effective adiabatic index, 7 m , such that 
B 2 oc p 7m . Possible values of this index are bounded by 7 m = for compres- 
sion/rarefaction parallel to the field lines of an ordered field and 7 m = 2 for 
compression/rarefaction perpendicular to the field, whereas 7 m < 1 seems to 
be indicated by observations. The order-of-magnitude increase in sound speed 
between the ionized and neutral gas means that the plasma /3-parameter (ra- 
tio of gas pressure to magnetic pressure) will be between 5 times (if 7 m = 0) 
and 2000 times (if 7, n = 2) higher in the Hll region than in the undisturbed 
neutral gas, assuming the Hll region to be five times overpressured. Thus, 
it is plausible that the magnetic field should play a much lesser role in the 
ionized gas than in the neutral gas. Nevertheless, the magnetic fields can still 
have dramatic effects on Hll region dynamics, particularly for the ionization 
front, and magnetohydrodynamic models of these are presented in detail in 
the chapter by Williams chapter in this volume. 



This is not the case for radiation pressure on the gas because of its dependence on 
the product of the flux and neutral density, which is roughly constant throughout 
the ionized region. 
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Instabilities 

The dynamics of the photoioni zed gas will als o be influenced by many dif- 
ferent kinds of instabilities fsee I Williams! |2003 for an overview), which may 
contribute to the detailed structure observed in H II regions. Different modes of 
instability at the ionization front are fo und, depending on whet her the shell of 
shocked neutral gas outside it is thick l)Williamsl E"002'l or thin l)Giuliani|ll97fll: 
iGarcia-Seeura fc Franco! Il996|) . Recom binations in the fully ionized gas are 
found t o damp the inst abilities in some l|AxforJl96J:lMizuta et al"1l2005|) but 
not all llWilliamsll20ni cases. The interaction between streams of gas inside 
the H II region, such as stellar winds or photoablation flows, provides oppor- 
tunities for further instabilities (e.g., Rayleigh- Taylor, Kelvin-Helmholtz) and 
these may interact with ionization front i nstabilities in complicat ed ways. An 
example is shown in the simulations of lArthur fc Hoard l)2005|) . where the 
fragmention of a wind-driven shell triggers shadowing instabilities of the type 
discussed hv lWilliamij lll999r> . 

3 H II region evolution 

While the classical scenario for the expansion of an ionized Stromgren sphere 
in a constant density medium has clear didactic value, it not necessarily rel- 
evant to the evolution of real Hll regions. Massive stars form in dense and 
highly inhomogeneous molecular clouds and begin to emit ionizing photons 
while they are still accreting mass, possibly via a disk (see the chapter by 
Franco & Hoare chapter in this volume for an overview of high-mass star for- 
mation). Various mechanisms may act to prolong the duration of the early 
phases of evolution, either by physical confinement or by providing a reservoir 
of neutral gas. Eventually the Hll region may break out from the molecular 
cloud and become optically visible, but its evolution will still be strongly af- 
fected by the strong density gradients in its environment. Since high-mass star 
formation is strongly clustered, extended Hll regions will tend to be excited 
by multiple OB stars. As the region's age begins to exceed the main-sequence 
lifetime of the highest mass stars, then the powerful winds from evolved stars 
and subsequent supernova explosions will have a dramatic effect on the dy- 
namics. 

3.1 Early phases: hypercompact and ultracompact regions 

The smallest observed Hll re gions are of size ~ 10~ 3 par sec, and are classi- 
fied as hypercompact regions l|Carri d et al .1119971: Scwil o et al.l At these 
small sizes, the escape velocity at the ionization front is larger than the ionized 
sound speed, and so gravitational effects are important. If the central star is 
still accreting via a spherical Bondi-Hoyle flow, then the accretion velocity 
just outside the ionization front exceeds the R-critical velocity of 2ci, so the 
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gas suffers onl y a mild de celeration in the front and continues its accretion 
onto the star l)Ketoll200^ . During this phase, the ionization front does not 
expand dynamically, but can only increa se its radius if the ionizing luminosity 
of the central star increases llKetoll2003l) . This phase is therefore rather long- 
lived, with a duration given by the timescale for the nascent star to grow by 
accretion, of order 10 5 years. 

Hll regions with sizes ~ 0.01-0.1 parsec and densities > 10 5 cm -3 are clas- 
sified as ultracompact and these are much more numerous and well-studied 
than the hypercompact r egions (e.g., IWood fc Chiirchwelllll98i iKurtz et alJ 
Il994t lde Pree et al.ll2005|) . They show a diverse range of morphologies, such 
as cometary, shell-like, bipolar, or irregular, although roughly half are spher- 
ical or unresolved. The high number of observed ultracompact Hll regions in 
the galaxy has been taken to imply a lifetime for this phase of ~ 10 6 years 
llWood fc Churchwelll H9891 , al though it should be note d that other estimates 
are as short as 5 x 10 4 years (jComeron fc Torralll996|) . A remarkable diver- 
sity of dynamical models have been proposed in order to explain this lifetime, 
which is longer than the time for a classical Stromgren sphere to expand past 
0.1 parsec. In some models, the ionized gas is physically confined to ultra- 
compact sizes, either by the ram pressure of the ambient medium as the star 
moves through the molecular gas llvan Buren et alJll99Cl>. or by an accretion 
flow onto the star (|Ketdl2002l iGonzalez-Aviles et al.l l2005). or merely by the 
high static pressure (thermal plus magnetic plus turbulent) o f the molecular 
cloud core I de Pree et al.lll99"itlGarcia-Segura fc Francolll99fil) . Other models 
allow the ionized gas to expand freely but provide a reservoir of neutral gas, 
the photoablation of which sustains the brightness o f the ultracompact core : 
either an accretion disk around the high-mass star ij Hollenb ach et alJ 1 19941 : 
lYorke fc Welzlll996t UTich ing fc Yorkdll997HLugo et alJl2004). or dense neu- 
tral clumps l|Dvson et alJll995l iLizano et al.lll996l iRedman et alJll99^ . In 
many instances, ultracompac t Hll regions are se en to be embedded in larger- 
scale, more diffuse emission ijKurtz et al.lll999|) . which would tend to favor 
the second class of models if the extended emission were physically connected 
with the ultracompact region. However, it is entirely possible that each of 
these models is applicable to some subclass of ultracompact regions. 



3.2 Later phases: compact and extended regions 

At least some ultracompact Hll regions eventually expand to larger sizes, 
forming the class of compact Hll regions l|Mezger et al.lll96^) . with sizes 0.1- 
0.5 parsec and densities ~ 10 4 cm~ 3 . Most compact regions are embedded 
inside extended Hll regions, which are typically several parsecs in size, with 
densities ~ 100 cm~ 3 . The di stribution of morphologi cal types is similar to 
that for ultracompact regions l(Hunterlll992t lFichlll99^ , although most mod- 
eling effort has gone into explaining the cometary shaped regions. These are 
similar in appearance to many op tically visible Hll regions, which are clas- 
sified as blister- type l|lsrae]lfl978l) . and of which the best-known example is 
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the Orion Nebula |o The low extinction to these optical regions 

proves that they must be on the near side of any accompanying molecular 
gas, and they generally show blueshifted velocities of order the ionized sound 
speed in optical emissio n lines, indicative of flow away from the molecular 
cloud l|Z ucker man! 1 1 9 73^1 . A broad class of models, commonly referred to as 
champagne flows, has been proposed to account for these objects. These mod- 
els share the property that strong density gradients in the neutral/molecular 
gas allow the ionization front to break out in some directions, leading to a 
flow of high-pressure ionized gas in the same direction. The original model of 
Tenorio- considered the one-dimensional propagation of an ion- 

ization front inside a dense cloud with a sharp edge. Once the ionization 
front reaches the edge of the cloud, it rapidly propagates through the much 
rarer intercloud medium and is followed by a strong shock that is driven by 
the higher pressure ionized cloud material. A rarefaction wave simultaneously 
travels back into the ionized cloud, initiating the champagne flow that ac- 
celerates the ionized gas up to several times the ionized sound speed as it 
flows away from the cloud . Following studies extended this work to two di- 
mensions (see references in lYorkejllQSd. section 3.3) and considered the effects 
of disk-shaped clouds and more gradual cloud boundaries. 

Early work on the dynamics of champagne flows considered scenarios that 
were intrinsically non-steady. Howe ver, it is also possible to construct quasi- 
stationary champagne flow models (jHennev in which the struc- 
ture of the ionized flow remains constant over several times its dynamical 
timescale. These models are valid so long as the density of the neutral gas 
that confines the ionization front on one side is high enough that the ioniza- 
tion front moves slowly and encounters constant upstream conditions during 
the evolution. It can be shown that in the steady-state limit a divergence of 
the ionized flow is necessary in order to produce significant acceleration. This 
is beacuse in a strictly plane-parallel geometry, and in the absence of body 
forces, the acceleration is proportional to the gradient of the sound speed, 
which is almost zero in the isothermal ionized gas. Examples of such flows 
with different degrees of divergence are shown in Fig. [21 The degree of diver- 
gence is controlled by the neutral density profile in the lateral direction (per- 
pendicular to the sharp density gradient at the surface of the cloud). When 
this density profile is constant, as in the upper panels of Fig. EI the ionization 
front is concave (negative curvature), leading to weak divergence of the flow 
and slow acceleration. For a steep lateral density profile, as in the lower panels 
of Fig. |21 the ionization front becomes convex (positive curvature), giving a 
strong, almost spherical, divergence to the flow, which now accelerates much 
more strongly. This configuration is similar to that seen in globule flows. Ap- 
proximate analytic calculations suggest that the ionization front should be 
flat when the lateral density profile is proportional to 1/(1 + r 2 ), in which 
the distance, r, from the symmetry axis has been scaled by the axial offset 
between the star and the ionization front. Numerical simulations confirm that 
this is indeed the shown in the central panels of Fig. |21 
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Fig. 2. Simulations of quasi-steady champagne flows, produced by the ioniz ation of 
the interface of a dense molecular cloud (adapted from lHennev et al.ll2005afl . Three 
different models (top to bottom) are shown with different lateral density profiles in 
the neutral gas. The left panels shows the ionized density, while the right panels show 
the total pressure with superimposed streamlines, the darkness of which indicates 
the gas velocity. The position of the ionizing star is shown by a black circle in each 
case. 



The stellar wind from the ionizing star may also be expected to have an 
effect on t he dy namics of the champagne flow and this was first modelled by 
IComeronl l|l997|) . An alternative explanat ion for the cometary sh apes shown 
by many regions is the bowshock model Ivan Buren et alj l|l99Cfl previously 
discussed in the context of ultracompact regions. In this model, the ionizing 
star moves supersonically through the molecular gas and the ionization front 
is trapped in the dense shocked shell formed by the interaction between the 
stellar wind and th e amb ient gas. An extension of this scenario was stud- 
ied bv lFranco et alJ l|20fl.4> . in which the ionizing star moves in and out of a 
high density molecular cloud core. An exhaustive study of the interplay be- 
tween champagne flow s, stellar winds, and stellar motion was carried out by 
lArthur fc Hoarel ([2005) , who found a great variation in the re sultant morphol- 
ogy, depending on the strength of the stellar wind. All the lArthur fc Hoarel 
models have concave ionization fronts, as in the top row of Fig. [3 If the stellar 
wind is weak, as in the first panel of Fig.0 then the champagne flow is hardly 
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Fig. 3. Simulations of the interaction of a champagne flow with a stellar wind (from 
lArthur fc HoareH2005Tl . Grayscale shows log of ionized density and arrows show gas 
velocity on two different scales: the longest white arrow represents 30 km s _1 while 
the longest black arrow represents 2000 km s _1 . The left panel shows a model with a 
weak wind (mass loss rate of 1O _7 M0 yr _1 ) and the right panel shows a model with 
a strong wind (1O _6 M0 yr _1 ). Both models have a wind velocity of 2000 km s _1 
and an ionizing photon luminosity of 2.2 x 10 48 s _1 . 



affected, except in a narrow region around the axis, where it forms a weak 
bowshock around the stellar wind. A stronger stellar wind, as in the second 
panel of Fig. [21 has a much more dramatic effect. The point of pressure bal- 
ance along the axis between the stellar wind and the H II region gas is pushed 
closer to the ionization front than where the sonic point would have been in 
the champagne flow. As a result, there is no champagne flow on the axis, al- 
though away from the axis a transonic champagne flow does develop, which 
then shocks against the stellar wind. The large obstacle provided by the wind 
means that the champagne flow acceleration is predominantly parallel to the 
ionization front in this case. The inclusion of a stellar velocity in the direction 
of increasing gradient produces only small changes in the results unless the 
stellar velocity is higher than the ionized sound speed. Such models always 
show a greater resemblance to champagne flow models than to classical bow- 
shock models unless the density gradient in the neutral cloud is very shallow. 
Either cometary or shell-like morphologies can be produced, depending on the 
steepness of the neutral density distribution. 

Another class of models considers the radial expansion of an H II region in 
a spherical cloud with a density distribut ion steeper than r~ 3 / 2 l|Franco et alJ 
ll99Ci IShu et~alll2002t iLizano et alJl200.^ . in which no static equilibrium so- 
lution exists for the position of the ionization front. When the star turns on, 
an R-type ionization front quickly propagates to large radii, and the ionized 
gas adopts a configuration in which a roughly uniform density core, bounded 
by a shock, expands slightly faster than the ionized sound speed. Such a re- 
gion is doomed to rapidly decline in brightness after a few sound crossing 
times. Confusingly, these entirely density-bounded models are also referred to 
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as champagne flows, although they are significantly different from the models 
considered in the previous paragraph. A classical champagne flow is ionization- 
bounded on at least one side and will maintain a high brightness indefinitely, 
as long as sufficient neutral gas exists to confine the ionization front close to 
the ionizing star on that side. 

At even larger scales, any expansion of the Hll region will stop when it 
comes into pressure balance with the local interstellar medium. The mean 
pressure in spiral ar ms at radiu s of solar circle is ~ 5 x 10 K cm -3 , which 
is 90% non-thermal ijCoxl 12005). This would balance that of the ionized gas 
when its density has fallen to a few cm~ 3 . The equivalent Stromgren radius 
is ~ 30 parsec for a typical O star, but would be would be larger for re- 
gions ionized by a cluster of stars, in which case vertical density gradients 
in the disk may become important. However, the time required for an Hll 
region to expand to such sizes is comparable with the evolutionary timescales 
for high-mass stars, so that the powerful winds from the later stages of stel- 
lar evolution and the expansion of supernova remnants will have a profound 
effect on the dynamics of the region. Simila r considerations apply to H II re- 
gions seen in external galaxies l|Shieldsll9 90). although t he largest H II regions 
(> 500 parsec) show highly supersonic velocities (e.g., iTerlevich fc Melnicld 
ll98lUGallagher fe Hunte^ll983tlRelafio et allEx)^ . which correlate with the 
luminosity of the region and may be indicative of virial equilibrium of the 
ionized gas in the gravitational potential of the star cluster. 



3.3 Clumps and turbulence 

One shortcoming of the models mentioned in the previous section is that 
they consider only smooth density distributions for the neutral gas into 
which the Hll region propagates. In reality, this is unlikely to be the case, 
since the molecular clouds in which high-mass stars form are known to 
be highly structured, wh ich is probably the result of supersonic turbulence 
(Elm egreen fc Scalol2004l) . Dynamic models that attempt to take into account 
the clumpy na ture of the neutral gas were first proposed by John Dyson and 
collab orators l|Dvson et alJll995l Twiiliams et alJll996t iRedman et al.lll99l 
1998) in order to explain the morphology and lifetimes of ultracompact H II 
regions. In these models dense neutral clumps are assumed to be distributed 
throughout the ionized volume and to act as mass-loading sites for a stel- 
lar wind from the ionizing star. The models are agnostic with respect to the 
physics of the mass injection process, which is not modeled directly, but which 
may be ablation by photoionization or hydrodynamic processes. Instead, the 
flow is treated statistically under the assumption of prompt and complete 
mixing of the ablated gas with the stellar wind. Similar models have also 
been presented that treat the ablation process in more detail and explicitly 
allow for the finite lifetime of the clumps ijLizano et al.ll99"5lArthur fc Lizanol 
Il997|) . Both sets of models result in a shell-like morphology for the H II region 
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(if the clumps are distributed uniformly), bounded by a recombination front 
through which the ionized wind flows. 




Fig. 4. Three-dimensional simula tion of the photoevaporation o f a clumpy self- 
gravitating globule (adapted from iKessel-Devnet fc Burkertll2003t) . Left: Phase of 
maximum compression, 1.6 x 10 5 years. Center: Equilibrium cometary phase, 2.8 x 
10 5 years. Right: Fragmentaion phase, 4.2 x 10 5 years. The grayscale measures the 
absolute value of the difference between the column density of the ionized and neutral 
gas (white corresponds to the ionization front). Box size is 0.85 parsec in the left 
panel and 1.8 parsec in the other panels. 

A more detailed treatment of the interaction of ionizing photons with a 
clumpy medium is possible only with three-dimensional numerical simula- 
tions. The first such simulation to be carried out studied the response of a 
static equilibrium self-gravitatin g spherical globule to an external ionizing flux 
ijKessel-Devnet fc Burkertl l2003h The density distribution of the initial glob- 
ule (of mass 40 M® and radius 1 parsec, on the margin of Jeans instability) 
is perturbed by Gaussian fluctuations of amplitude ~ 50%. It is found that 
the turbulence generated behind the shock driven into the globule by the ion- 
ization front is sufficient to prevent gravitational collapse of the globule at its 
point of maximum compression, despite collapse having occurred at this point 
in a companion simulation of a globule without density perturbations. After 
~ 5 x 10 5 years, the neutral gas that has survived photoablation has broken 
up into smaller sub-globules (Fig.^J. Similiar simulations have also been car- 
ried out of the photoevaporatio n of much smaller, non-self-gravitating, clumpy 
globules l|Gonzalez et al.ll2005^) . In this case, although some fragmentation of 
the globule occurs, the fragments remain physically attached to one another 
until the globule is completely evaporated. 

Simulations of hydro dynamical and magnetohydrodynamical turbu lence 
in molecular clouds fe.g.. lMac Lowlll99^ : IVazauez-Semadeni et a l~ll20051 have 
now reached a level of refinement that makes them su i table to use as ini- 
tial conditions for the evolution of Hll regions. iLi et al.l (|Jo04) carried out a 
preliminary investigation of this problem by studying only the initial R-type 
propagation of the ionization front, before the dynamics of the ionized gas 
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Fig. 5. Dyn amical evolution of a n Hll region in a turbulent molecular cloud, 
adapted from IMellema et al.l <l2005f) . For each of two evolutionary times, 100,000 
and 300,000 years, five yz slices through the (4 parsec) 3 computational box are 
shown, with offsets from the ionizing star as shown on the bottom axis. Density 
is shown on a negative logarithmic scale between 10 (white) and 10 5 cm 3 (black). 
Temperature is shown on a positive linear scale between (black) and 10 4 K (white) . 



becomes important. This approach was extended bv IMellema et all l)2005|h 
who carried out a full radiation-hydrodynamic simulation of the evolution of 
an Hll region in a turbulent medium, following the birth of an ~ 25 Mq 
star inside the densest molecu lar clump (~ 10 6 cm -3 ) that formed in the 
IVazauez-Semadeni et al turbulence simulation. 

The results of this simulation are shown in Fig. [3] for two evolutionary 
times. At the earlier time, the Hll region is mainly confined to a compact 
core of radius ~ 1 parsec, although the ionization front has already broken 
out to the boundary of the grid through a corridor of low-density neutral gas 
in one direction, as can be seen in the lower-right corner of the rightmost 
temperature image. The density variations in the ionized gas are rather mild 
compared with those in the neutral gas and are of a different nature. Instead 
of dense sheets and filaments, one finds an almost constant density gas filling 
roughly half of the ionized region, with the remainder occupied by low-density 
cavities c arved out by transonic photoablation flows, as described in section l2~H 
(see also IWilliams et al"1l200ll lHennevl l2003:). This structure of the ionized 
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gas can be better appreciated at the later evolutionary time, where multiple 
photoablation flows are seen, streaming off dense neutral condensations. The 
highest ionized densities are found at the ionization fronts of these flows and 
their mutual collisions produce many weak shocks, which generate density 
structure in the interior of the Hll region. Some of these structures become 
compressed sufficiently to significantly absorb the ionizing photons, causing 
the ionized gas beyond them to temporarily recombine. Regions where this 
has occured can be seen as intermediate-temperature gas in the bottom row 
of Fig. [SJ The photoablation flows help maintain a high velocity dispersion 
of the ionized gas, which remains roughly equal to the ionized sound speed 
during the entire evol ution, even though th e net radial expansion velocity falls 
to much lower values (Mcllcma ct al. 2005, Fig. 4). This mechanism is similar 
to that proposed bv lDvsorJ l| 19681) to explain the velocity dispersion observed 
in the Orion Nebula. 




Fig. 6. Fake emission line images from the simulation shown in Fig. |E|at an evolu- 
tionary time of 400,000 years. Top row: [Nn] 6584 A. Bottom, row: [Om] 5007 A. 
Views along the x, y, and z axes are shown from left to right, respectively. 

The predicted appearance of these simulations in optical emission lines 
is shown in Fig. H3 The [N n] line traces the low-ionization gas close to the 
ionization front, and these images are dominated by the bright rims of pho- 
toablation flows, while the [Om] line traces the higher ionization gas in the 
interior of the nebula, which shows a more diffuse aspect in the images. The 
predicted emission line structure is very similar to that observed in real neb- 
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Fig. 7. The Lagoon Nebula (M8). Combination of red and infrared continuum im- 
ages from the Second Epoch Sky Survey of the UK Schmidt Telescope, operated by 
the Anglo-Australian Observatory, digitized by the Space Telescope Science Insti- 
tute. A highly non-linear grayscale mapping has been used to show both bright and 
faint structures at the same time. The field of view is ~ 27 x 21 parsec, assuming a 
distance of 1500 parsec. 

ulae, such as the Lagoon Nebula shown in Fig.0 The broad-band niters used 
in that image include optical emission lines of both low and high ionization, 
so it should b e compared w ith the superposition of the [N n] and [O m] images 
(see Fig. 1 of lDufourlll99i for pure emission line images of the inner regions 
of this nebula). As in the simulations, one sees sharply bounded bright arcs 
and dark extinction features in the low-ionization periphery of the nebula, 
together with more diffuse emission in the higher ionization core. It should be 
noted that the Lagoon Nebula is several times larger than the simulation and 
is ionized by several massive stars, but a similar morphology is seen in more 
compact nebulae, such as M16, M20, and M42. 

Such sculpted structures at the boundary of the nebula can also be pro- 
duced via instabi lities, even if the ambient neu tral gas is initially smoothly 
distributed (e.g., iGarcia-Segura fc Francdll99r^ . A clear example of this is 
seen in Fig. 3 of IWilliamsl (2003), where large-scale density gradients lead 
to a thin-shell instability and the formation of multiple photoablation flows 
from the now clumpy shocked neutral layer. Whether ionization front insta- 
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bilities or pre-exisiting structures in the neutral gas are more important in a 
given region is still an open question. However, the density contrasts found 
in molecular clouds on parsec scales are much more extreme than those seen 
in Hll regions, due to the highly supersonic nature of the turbulence and 
the effects of self-gravity. Thus, it seems likely that the structure of the neu- 
tral gas is dominant in shaping Hll regions at these scales. At smaller scales, 
the turbulent velocity dispersion, and hence the density contrast, is much re- 
duced, so that ionization front instabilities may play a more important role, 
although these tend t o saturate at amplitudes of order the recombination 
length 1williamsl Eo02h which is approximately ten times the ionization front 
thickness, or ~ n~ x parsec, where n is measured in cm~ 3 . 

Summary 

Hll regions, particularly the extended, optically visible variety have the repu- 
tation of being "messy" objects, ill-suited to the austere theoretical prejudice 
of perfect symmetry. However, as I have tried to illustrate in this review, even 
their messiest aspects are beginning to fall under the purview of dynamic 
modelling. We are still some way from a totally satisfactory model of any par- 
ticular object, let alone Hll regions as a class, but a large amount of progress 
has been made. 
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